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^\. . A quantitative model describing large magnetostrain effect observed in several ferromagnetic 

^vj ' shape memory alloys such as Ni2MnGa is briefly reported. The paper contains an exact thermody- 

namic consideration of the mechanical and magnetic properties of a similar type materials. As a 
result, the basic mechanical state equation including magnetic field effect is directly derived from 
^ ' a general Poisson's rule. It is shown that the magnetic field induced deformation effect is directly 

^ , connected with the strain dependence of magnetisation. A simple model of magnetisation and its 

dependence on the strain is considered and applied to explain the results of experimental study of 
large magnetostrain effects in Ni2MnGa. 
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In addition to some giant magnetostriction materials, ferromagnetic shape memory alloys was recently suggested 
as a general way for the development of a new class of the magnetic- field-controlled actuator materials ||l|,||,|[^,|l^,0 . 
It is now a goal of research projects in several groups directed on the development of ferromagnetic alloys exhibiting 
fj ■ also a martensitic phase transition that would allow control of large strain effect by application of a magnetic field at 
constant temperature in a martensitic state. Numerous candidate shape memopy materials were explored including 
Ni2MnGa, Co2MnGa, FePt CoNi, and FeNiCoTi during past few years iq,H,nl. Magnetically driven strain effect is 
expected to occur in these systems. According to results reported in g] the large strains of 0.19% can be achieved in 
/y-\ i a magnetic field of order 8 kOe in the tetragonal martensitic phase at 265 K of single-crystal samples of Ni2MnGa. 
fT^ ' This strain is an order of magnitude grater than the magnetostriction effect of the parent, room temperature cubic 
■TJ" , phase. 

^^ ' Ni2MnGa is an ordered by L21 ferromagnetic Heusler alloy having at high temperature cubic (a = 5.822 A) 

^P [ crystal structure that undergoes martensitic transformation at 276 K into a tetragonally distorted structure with 

crystalline lattice parameters: a=b=5.90 A and c=5.44 A |6|. The martensitic phase accommodates the lattice 

distortion connected with transformation by formation of three twin variants twinned usually on (110) planes and 

jrt ' having the orientation of the tetragonal symmetry axes nearly to three possible [100] directions. The saturation 

value of magnetization was found to be about of 475 G. Magnetization curve of the low-temperature twinned phase 

usually displays two-stage structure at 265 K ||] with a sharp crossover at about of 1.7 kOe from easy low-field 

'■^ ■ magnetization below to a hard stage above this value up to the 8 kOe saturation field value. Such a behavior is 

^ [ connected with a different response of different twin variants to the applied field. The measurements usually show 

O ■ a definite magnetostrain value along [100] as a function of the magnetic field applied in the same direction [g|. It is 

generally expected that a large macroscopic mechanical strain induced by the magnetic field in similar type systems 

is microscopically realized trough the twin boundaries motion and redistribution of different twin variants fractions 

in a magnetic field. The main thermodynamic driving forces have in this case magnetic nature and connected with 

, ; high magnetization anisotropy and differences in magnetization free energies for different twin variants of martensite 

The main goal of this brief publication is to give the right thermodynamic consideration of the mechanical and 
magnetic properties of a similar type materials and represent the quantitative model describing large magnetostrain 
effect observed in several ferroelastic shape memory alloys such as Ni2MnGa. It is shown that the magnetic field 
induced deformation effect directly follows from the general thermodynamic rules such as Poisson equation and 
connected with the strain dependence of magnetization. A simple model of magnetization for the internally twinned 
martensitic state and its dependence on the strain is considered and applied to explain the results of experimental 
study of large magnetostrictive effects in Ni2MnGa. 
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II. GENERAL THERMODYNAMIC CONSIDERATION 

Consider the general thermodynamic properties of the materials which can show both the ferroelastic and the 
ferromagnetic properties. Most of the shape memory alloys usually display ferroelastic behavior in martensitic state 
connected with redistribution of different twin variant fractions of martensite under the external stress applied through 
the motion of twin boundaries. Ferromagnetic shape memory materials have an additional possibility to activate 
the deformation process in a twinned martensitic state by the application of magnetic field simultaneously with 
magnetization of the material. According to general thermodynamic principles both the mechanical and the magnetic 
properties of similar type materials can be represented by the corresponding state equations: 

a — a (e, /i) (1) 

m ~ m{e, h) (2) 

where, the Eq.™) reflects the mechanical properties through stress-strain a — e equation in presence of magnetic field 
h and Eq.(p[) gives the magnetization value m as a function of magnetic field applied h and strain e. Both these 
equation can be obtained from an appropriate thermodynamic potential as follows: 

a{e,h)^—G{e,h) m{e,h) ^ -—G{e,h) (3) 

where, G (e, h) = G (e, h) — hm{e, h), and G (e, h) is the specific Gibbs free energy at fixed temperature and pressure 
condition. Both state equation are not completely independent functions and must satisfy known Poisson's rule: 

-aie,h) = --mie,h) (4) 

Integration of this equation over the magnetic field starting from /i = at a fixed strain gives an important represen- 
tation of the mechanical state equation including magnetic field effects: 

h 

d f 

a — ao (e) — -— / dhm (e, h) (5) 

oe J 



According to this equation the external stress on the left is balanced in equilibrium by both the pure mechanical stress 
ctq (e) = cr (e, 0) resulting from the mechanical deformation of the material ai h = and the additional magnetic field 
induced stress that is represented by the second term on the right in this equation. It is also important to note 
that all the effect of the magnetic field on the mechanical properties is directly determined by the strain dependence 
of magnetization. In a particularly important case: a = const = one can obtain a general equation determining 
magnetically induced strain (usually called as a magnetic shape memory or MSM-effect) as follows: 



n 

d f 

CTo (e) = — / dhm (e, h) 



(6) 



and its linearized solution: 



de 



e=0 




W=(^) \^ I dhm{e,h)\ (7) 



that can be used when e is much less than a martensite lattice tetragonal distortion value eg- According to Eqns.(y) 
and (m) the magnetization and its dependence on the strain is responsible for the MSM-effect and is the main subject 
for detailed discussion and modeling. 



III. THE MODEL AND ITS APPLICATION TO NI2MNGA TETRAGONAL MARTENSITE 

Consider a typical situation corresponding to measurements of large strain induced by the magnetic field in the 
tetragonal internally twinned martensite of Ni2MnGa obtained from the austenitic single crystal studied in |^ when 
the magnetic field is applied along [100] direction of parent austenitic phase and strain measurements were performed 
in the same axial direction. In this case the crystallographic [100] [010] [001] axes for all three possible twin variants 
of the tetragonal martensitic phase will be nearly parallel to the external field applied. More exactly, the crystallo- 
graphic orientation relationships between the austenitic and martensitic phases can be obtained by using the usual 
methods of the crystallographic theory. Additional small rotations of the tetragonal phase axes are expected but the 
corresponding rotation angles can not exceed few degrees in the case of Ni2MnGa and may be neglected for simplicity. 
Fig.l. schematically shows the expected alignment of the magnetic field applied, crystallographic orientations and 
magnetization curves for three possible tetragonal phase variants. 
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FIG. 1. Schematic representation of the magnetic field alignment and magnetisation behavior for three different crystallo- 
graphic variants of the tetragonal martensitic phase. 

Therefore, the magnetic field is applied along the tetragonal symmetry axis only for one type twin variant (which is 
called here the axial, or a-type) and simultaneously in the transversal direction in respect to the tetragonal symmetry 
axes of another two (transversal, or t-type) variants. The investigation of magnetization properties of Ni2MnGa 
performed for a single tetragonal variant of martensite obtained by the mechanical compression method |12] has 
shown a considerable difference between the magnetization curves along the tetragonal [100] direction in comparison to 
another transversal [010] and [001] directions. It was found that a tetragonal axis is the easiest magnetization direction 
and requires considerably less value of saturation field ha than a saturation field ht characterizing magnetization in 
two hard transversal directions as it is schematically shown in Fig. 1 . In a general case the calculation of magnetization 
for the material with a complicated twin microstructure geometry requires some special approach. In this paper we 
will ignore, for simplicity, a similar type problem and will consider these effects in other publications. Taking into 
account the presence of magnetic anisotropy and difference in magnetization behavior between the axial ma{h) and 
transversal mt{h) twin variants we consider a simple model of magnetization for the multi- variant martensitic state 
that gives the main contribution into magnetisation insensitive to the fine details of twin microstructure. This model 
treats the multi-twin martensitic state as a composite material consisting of an easy magnetization area occupied 
by axial type twins and hard magnetization region of two transversal twin variants. Denoting as x the total volume 
fraction of the axial twin domain and (1 — x) - transversal type twin domain fractions, respectively, one can write the 
magnetization of the material as follows: 

m {x, h) — xma{h) + (1 — x)mt{h) (8) 

where, ma{h) and mt{h) are specific magnetization functions for the axial and transversal variants, respectively. On 
the other hand, the macroscopic strain along the axial direction can be found from a similar type equation: 

e = xel + {l-x)e",^^eo{x~^) (9) 

where, the diagonal matrix elements e" — £q and e^ — —3^0 represent the relative tetragonal distortion of the 
martensite crystal lattice along its tetragonal axis and two transversal directions, respectively. The compression 
distortion eo = 5.4% along the tetragonal symmetry axis was found in case of Ni2MnGa. One can easily exclude the 
fraction dependence from these two equations and obtain the magnetization as function of the macroscopic strain for 
the internally twinned martensitic state: 

(1 2 12 

TO (e, h) = <^ -^maih) + -mt{h) ^ + - (e/eo) {maih) - rntih)} (10) 

This equation immediately reproduces all the main peculiarities of the experimental magnetization curve including 
the sharp change of its slope at h = 1.75 kOe as indicated in Fig. 2. This singularity exactly appears at h = ha 
where the easy stage of magnetization process inside of the axial twin variants domain is finished. According to 
Eq.(|lO|) TO (£0, h) = rua (h) and to (— eo/2, h) = rrit (h) so, one can use this fact to obtain both the ha = 1-75 kOe and 
h — 8 kOe from the experimental magnetization curves measured in the multi- variant state. The model magnetization 
curve to (0, ft.) corresponding to zero strain value shows the same type behavior and singularity in slope as the 



experimental one. The difference between them is caused by the second term in Eq.(|10|) that gives an additional 
strain dependent contribution into the magnetisation. This contribution is directly connected with the MSM-effect 
and can be easily taken into account just after its calculation. 
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FIG. 2. Experimental and model magnetisation curves that show the compression and tensile strain effect on the macroscopic 
magnetisation of martensitic phase caused by the relative change of the twin variant fractions. 

One can also obtain the final equations representing the effect of magnetic field on the strain by using basic 
expressions ([7|) derived before from the general thermodynamic consideration: 



-1 '* 
(/i) = - f eo-p) dh {ma{h) ~ mt{h)} 



(11) 



IV. DISCUSSION AND CONCLUSIONS 

As follows from this equation, two factors determine the strain value and its field dependence. The first one 
is proportional to the slope of stress-strain curve and can be found from the usual mechanical compression test 
without magnetic field applied. The integral term reflects the effects of magnetic anisotropy and determines the 
functional magnetic field dependence of the strain. In particular, in absence of the magnetization anisotropy when 
ma{h) = mt{h) deformation effect is also vanishes. The saturation level of the strain is achieved at h = ht and above 
where ma{h) = mt{h) = mf^* and where the material has its maximal magnetization value ?7i*°*. One can easily 
obtain the corresponding saturation value of the strain performing the necessary integrations in Eq. (|l l[) as follows: 

e7:r = \ (eo^) _ {ht - K) m-* (12) 

Precise quantitative estimation of the saturation strain requires, in general, the corresponding mechanical testing. 
Here, we will use a simple estimation of da^/de ^ cto/eo- So, e^t™ ~ \ {<^o)~ (^t — ^a) ™'*''* where the characteristic 
stress (To representing ferroelastic mechanical behavior of the material is expected to be about of 20MPa in Ni2MnGa 
martensite. Using also the values of ht ~ 8fcOe , ha ^ l.75kOe and m''"* ~ 475G found from the magnetization curve 
analysis one can obtain a simple estimation: e^f^ ~ 0.49%. More precise estimation that follows from the mechanical 
testing results gives dao/de ~ (2-^3) ctq/so. Consequently, e™ j™ ~ (0.24 -;- 0.16)% which is in a better quantitative 
agreement with e™^™ ~ 0.14% experimental value. In order to achieve the larger magnetostrain effect comparable 
with the lattice tetragonal distortion value Eq ~ 5% one will need materials with very low ctq ^ 2MPa detwinning 
stress value. This task seems can be considered as the realistic one because the observation of cto ^ 2MPa [g| and 
(To ~ 8MPa ||l^] were reported in some publications. 

Fig. 3. shows the field behavior of the strain that follows from the model and its change for the different values of 
the magnetic anisotropy factor k — ha/ht defined as a ratio between the axial and transversal saturation fields. 
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FIG. 3. Magnetic anisotropy effect on the strain vs. magnetic field behavior according to the model calculations. 

The dimensionless strain response e™*™ (/j) /emax normalized by, 

ema. = I (eo^) ' /itm-* (13) 

increases from zero value at fc = 1 simultaneously with a corresponding shape change and shows the maximal possible 
deformation effect and linear type singularity for the low field strain behavior ai k — 0. This case corresponds 
to the maximally strong anisotropy when the axial saturation field becomes infinitely small ha — > and ma{h) 
immediately achieves its saturation level tti'"'* starting from an arbitrary low magnetic field and then still remains 
equal to a constant m*"* value during the magnetization process. Therefore, one can conclude that a linear low field 



behavior usually predicted in some previously developed models |l^Jll|l is directly connected with their assumption 
on the complete saturation of magnetization for the axial type twin variants. According to the present model such 
a type of assumption can be physically reasonable in the limit ha — > only. In other case the strain shows the 
normal parabolic type behavior in the low field h < ha region in agreement with the experimental observations. A 
good correspondence between the model and experimental results is indicated in Fig. 4. We neglected here the small 
hysteresis effects which are usually observed assuming to give the more detailed discussion of this problem by using 
some new developments and quantitative descriptions of hysteresis in shape memory materials. 
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FIG. 4. Magnetostrain effect: comparison results between the model and the experiment 
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